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Abstract 

Ratio asymptotics for matrix orthogonal polynomials with recurrence coefficients 
A n and B n having limits A and B respectively (the matrix Nevai class) were ob- 
tained by Duran. In the present paper we obtain an alternative description of the 
limiting ratio. We generalize it to recurrence coefficients which are asymptotically 
periodic with higher periodicity, or which are slowly varying in function of a pa- 
rameter. Under such assumptions, we also find the limiting zero distribution of 
the matrix orthogonal polynomials, generalizing results by Duran-Lopez-Saff and 
Dctte-Reuther to the non-Hermitian case. Our proofs are based on 'normal family' 
arguments and on the solution to a quadratic eigenvalue problem. As an appli- 
cation of our results we obtain new explicit formulas for the spectral measures of 
the matrix Chebyshev polynomials of the first and second kind, and we derive the 
asymptotic eigenvalue distribution for a class of random band matrices generalizing 
the tridiagonal matrices introduced by Dumitriu-Edelman. 

Keywords: matrix orthogonal polynomial, (block) Jacobi matrix, recurrence 
coefficient, (locally) block Toeplitz matrix, ratio asymptotics, limiting zero distribu- 
tion, quadratic eigenvalue problem, normal family, matrix Chebyshev polynomial, 
random band matrix. 



1 Introduction 

Let (P n (x))™ =0 be a sequence of matrix- valued polynomials of size r x r (r > 1) generated 
by the recurrence relation 

xP n (x) = A n+1 P n+1 (x) + B n P n (x) + A* n P n ^(x), n>0, (1.1) 

with initial conditions P-\{x) = G C rxr and Pq(x) = I r , with I r the identity matrix of 
size r. The coefficients A^ and are complex matrices of size r x r. We assume that 
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each matrix is nonsingular and B^ 
Hermitian conjugation. 



is Hermitian. The star superscript denotes the 



The polynomials generated by (1.1) satisfy orthogonality relations with respect to 



a matrix- valued measure (spectral measure) on the real line (Favard's theorem [5]) and 
they are therefore called matrix orthogonal polynomials. The study of such polynomials 
goes back at least to [18J and we refer to the survey paper [SJ for a detailed discussion 
of the available literature and for many more references. Some recent developments and 
applications of matrix orthogonal polynomials can be found in [31 [9J, [T7] among many 
others. 



To the recurrence (1.1) we associate the block Jacobi matrix 



{B A 1 

A*, By An 



4* 



A 



n-l 



o \ 



An-l 
B n -i) 



(1.2) 



which is a Hermitian, block tridiagonal matrix. It is well-known that 



det P n {x) = c n det(xl n - J n 



with c n = det(A n 1 • • • A 2 1 A 1 ) ^ 0, see [8[ I14j. and with the zeros of the matrix 
polynomial P n (x) we mean the zeros of the determinant detP n (x), or equivalently the 



eigenvalues of the matrix J n defined in (1.2) (counting multiplicities). 

The polynomials (P n (x))^ =0 are said to belong to the matrix Nevai class if the limits 



lim A n 

n— >oo 



A, 



lim B n 

n— >oo 



B. 



(1.3) 



exist, where we will assume throughout this paper that A is nonsingular. 

One of the famous classical results on orthogonal polynomials is Rakhmanov's theo- 
rem [23], see [UEIET] for a survey of the recent advances in this direction. Rakhmanov's 
theorem for matrix orthogonal polynomials on the real line is discussed in [9[ [29] . These 
results give a sufficient condition on the spectral measure of the matrix orthogonal poly- 
nomials, in order to have recurrence coefficients in the matrix Nevai class, with limiting 
values A = I r and B = 0. 



Duran [13J shows that in the matrix Nevai class (1.3), the limiting matrix ratio 



R(x) :-- 



lim P n (x)P~^(x), xGC\ [-M, M], 



(1.4) 



exists and depends analytically on x £ C \ [-M, M\. Here M > is a constant such that 
all the zeros of all the matrix polynomials P n (x) are in [—M,M]. Moreover, it is also 
proved in [13] that R(x)A~ 1 is the Stieltjes transform of the spectral measure for the 
matrix Chebyshev polynomials of the second kind generated by the constant recurrence 
coefficients A and B; see Section [4j 
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The present paper has several purposes. First we will give a different formulation of 
Duran's result on ratio asymptotics [13]. In particular we will give a self-contained proof 
of the existence of the limiting ratio R(x) and express it in terms of a quadratic eigenvalue 
problem. Second, our approach can also be used to obtain ratio asymptotics for some 
extensions to the matrix Nevai class. More precisely, we will allow the coefficients A n 
and B n to be slowly varying in function of a parameter, or to be asymptotically periodic 
with higher periodicity. In both cases we prove the existence of the limiting ratio (the 
limit may be local or periodic) and give an explicit formula for it. 

To prove these results we will work with normal families in the sense of Montel's 
theorem. Similar arguments can be found at various places in the literature and our 
approach will be based in particular on the work by Kuijlaars-Van Assche [21] and its 
further developments in [31 El [HI ES] . We point out that an alternative approach for 
obtaining ratio asymptotics, is to use the generalized Poincare theorem (see |23} [26] for 
this theorem); but the normal family argument has the advantage that it also works for 



slowly varying recurrence coefficients; see Section 3.1 for more details 



The third purpose of the paper is to obtain a new description of the limiting zero distri- 
bution of the matrix polynomials (P n (x))^ =0 . Duran-Lopez-Saff [15J showed that in the 



matrix Nevai class (1.3), and assuming that A is Hermitian, then the zero distribution 
of P n (x) has a limit for n — > oo in the sense of the weak convergence of measures. They 
expressed the limiting zero distribution of P n {x) in terms of the spectral measure for 
the matrix Chebyshev polynomials of the first kind, see Section [5] for the details. In 
contrast to the work of [15] the results derived in the present paper are also applicable 
in the case when the matrix A is nonsingular but not necessarily Hermitian. Maybe not 
surprisingly, the limiting eigenvalue distribution of the matrix J n in (1.2)-(1.3) is the 
same as the limiting eigenvalue distribution for n 



oo of the block Toeplitz matrix 



A* 



A 
B 

A* 



A 



0\ 



••• A 

A* Bj 



(1.5) 



The eigenvalue counting measure of the matrix T n has a weak limit for n — > oo [28] and 
a description of the limiting measure can be obtained from j5[ [TU1 128| . 

In this paper we establish the limiting zero distribution of the matrix polynomials 
P n {x) as a consequence of our results on ratio asymptotics. We also find the limiting 
zero distribution for the previously mentioned extensions to the matrix Nevai class. That 
is, the coefficients A n and B n are again allowed to be slowly varying in function of a 
parameter, or to be asymptotically periodic with higher periodicity. 

Incidentally, we mention that it is possible to devise an alternative, linear algebra 
theoretic proof for the fact that the matrices J n and T n have the same weak limiting 



eigenvalue distribution, using the fact that the block Jacobi matrix (1.2) is Hermitian 
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|20| ; however our approach has the advantage that it can be also used in the non- 
Hermitian case, at least in principle. This means that most of the methodology derived 
in this paper is also applicable in the case when the recurrence matrix (1.2) generating 
the polynomials P n (x), is no longer Hermitian, or when it has a larger band width in its 
block lower triangular part, as in [3]. In fact the key places where we use the Hermiticity 
are in the proof of Proposition 2.1 and in the proof of Lemma 7.1 see |llj ; but it is 
reasonable to expect that both facts remain true for some specific non-Hermitian cases 
as well. This may be an interesting topic for further research. 



The remaining part of this paper is organized as follows. In the next section we state our 
results for the case of the matrix Nevai class. In Section [3] we generalize these findings, 
to the context of recurrence coefficients with slowly varying or asymptotically periodic 
behavior. In Section [4] we apply our results to find new formulas for the spectral mea- 
sures of the matrix Chebyshev polynomials of the first and second kind. In Section [5] we 
relate our formula for the limiting zero distribution to the one of Duran-Lopez-Saff |15| . 
In Section [6] we indicate some potential applications in the context of random matrices. 
In particular we derive the limiting eigenvalue distribution for a class of random band 
matrices, which generalize the random tridiagonal representations of the /3-ensembles, 
which were introduced in [12] . Finally, Section [7] contains the proofs of the main results. 



2 Statement of results: the matrix Nevai class 
2.1 The quadratic eigenvalue problem 



Throughout this section we work in the matrix Nevai class (1.2)-(1.3). Observe that the 
limiting ratio R(x) in (1.4) satisfies the matrix relation 



A*R(x) + B-xI r + AR-\x) = 0, x G C \ [-M, M] , 



(2.1) 



which is an easy consequence of the three term recurrence ( 1.1 ). For a simple motivation 
of our approach we assume at the moment that for each x G C \ [— M, M], the matrix 
R(x) is diagonalizable with distinct eigenvalues z k = z k (x), k = 1, . . . ,r. Let v k (x) G C r 
be the corresponding eigenvectors so that 

R(x)\ k (x) = z k (x)v k (x), 



for k = 1, . . . , r. Multiplying (2.1 ) on the right with \ k (x) we find the relation 

(z k (x)A* +B-xI r + Zk ~ \x)A) v k (x) = 0, (2.2) 

where denotes a column vector with all its entries equal to zero. This relation implies 
in particular that 

det (z k (x)A* +B-xI r + z~ 1 (x)A) = 0. (2.3) 



We can view (2.3) as a quadratic eigenvalue problem in the variable z = z k (x), and 
it gives us an algebraic equation for the eigenvalues of the limiting ratio R(x). The 
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corresponding eigenvectors Vk(x) can then be found from (2.2). Note that the equation 



(2.3) has 2r roots z = zj-{x), k = 1, . . . , 2r. Ordering these roots by increasing modulus 

< \zi(x)\ < \Z 2 {X)\ <...< \z r (x)\ < \Zr+l(x)\ < . . . < \z 2r (x)\, (2.4) 

then we will see below that the eigenvalues of R(x) are precisely the r smallest modulus 
roots zi(x), . . . , z r {x). 

In order to treat the general case of eigenvalues with multiplicity larger than 1 we 
now proceed in a more formal way. Inspired by the above discussion, we define the 
algebraic equation 

= f(z, x) := det(zA* + B + z~ x A - xl r ), (2.5) 



where z and x denote two complex variables. As mentioned, one may consider ( |2.5[ ) to 
be a (usual) eigenvalue problem in the variable x and a quadratic eigenvalue problem in 
the variable z. 



Expanding the determinant in (2.5), we can write it as a Laurent polynomial in z: 



f(z,x)= h{x)z\ (2.6) 

k=—r 

where the coefficients fk(x) are polynomials in x of degree at most r, and with the 
outermost coefficients f r (x) and f- r {x) given by 

f r (x) = f r = det A* , /_ r (a:) = /_ r = det A 

Solving the algebraic equation f(z,x) = for z yields 2r roots (counting multiplici- 
ties) Zk = Zk(x) which we order by increasing modulus as in ( |2.4[ ). If x G C is such that 
two or more roots Zf t (x) have the same modulus then we may arbitrarily label them such 



that (2.4) holds 



Define 

r := {x e C | \z r (x)\ = |«r + l(a:)|}. (2.7) 
It turns out that the set Tq attracts the eigenvalues of the block Toeplitz matrix T n in 



(1.5) for n — > oo, see [5"j [2"H]. This set will also attract the eigenvalues of the matrix 



J n in (1.2). The structure of Tq is given in the next proposition, which is proved in 



Section 17.11 

Proposition 2.1. Tq is a subset of the real line. It is compact and can be written as 
the disjoint union of at most r intervals. 

2.2 Ratio asymptotics 



For any x E C and any root z = Zk{x) of the quadratic eigenvalue problem (2.5), we can 



choose a column vector v^(x) G C r in the right null space of the matrix Zk{x)A* + B 



xl r + z k 1 (x)A, see fl!E2] ). If the null space is one-dimensional then we can uniquely fix 
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the vector v k (x) by requiring it to have unit norm and a first nonzero component which 
is positive. 

If the null space is 2- or more dimensional, then the vector v k (x) is not uniquely 
determined from (2.2). We need some terminology. (The following paragraphs are a 
bit more technical and the reader may wish to move directly to Theorem 2.3 ) For any 
x,z £C define d(z,x) as the geometric dimension of the null space in (2.2), i.e., 



d(z, x) := dim{v E C r | (A*z + B + Az' 1 - xl r )v = 0}. 
Also define the algebraic multiplicities 



mi(z, x) 
rri2(z, x) 



max{A: E Z> | (Z - z) k divides Z r f(Z,x)}, 



max{k E Z>o | {X 



divides f(z, X)}, 



(2.8) 



(2.9) 
(2.10) 



where Z and X are auxiliary variables and the division is understood in the ring of 
polynomials in Z and X respectively. 

In the spirit of linear algebra, we can think of d(z, x) as the geometric multiplicity 
of (z,x) E C 2 while m\(z,x) and rri2(z,x) are the algebraic multiplicities of (z,x) with 
respect to the variables z and x respectively. 

The next lemma is a generalization of a result of Duran [13., Lemma 2.2]. 

Lemma 2.2. (Algebraic and geometric multiplicities:) Let A,B E C rxr be matrices 
with A non-singular and B Hermitian. Define f(z,x) by (2.5). Then the algebraic and 
geometric multiplicities in (2.8)-( 2T0| ) are the same: 

d(z, x) = mi(z, x) = 7Tl2(z, x), 

for all but finitely many z,x E C. 



Lemma |2.2| is proved in Section [7. 2[ We note that the particular Hermitian structure 
of the problem is needed in the proof. 

Now let x E C and consider the roots zi(x), . . . , Z2 r (x) in (2.4) with algebraic multi- 
plicities taken into account. Lemma |2 . 2 1 ensures that for all but finitely many x E C, we 



can find corresponding vectors vi(x), . . . , V2 r (x) having unit norm, satisfying (2.2), and 
such that 

dim{v;(x) | / E {1, . . . ,2r} with z t (x) = z k (x)} = d(z k (x),x) = mi(z k (x),x), (2.11) 

for each fixed k. In what follows we will always assume that the vectors v k (x) are 
chosen in this way. We write ScC for the set of those x £ C for which (2.11) cannot 
be achieved. Thus the set S has a finite cardinality. 

Now we are ready to describe the ratio asymptotics for the matrix Nevai class. The 
next result should be compared to the one of Duran [13] • 
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Theorem 2.3. (Ratio asymptotics.) Let A,L>£ c r><r 5 e matrices with A non-singular 
and B Hermitian. Let P n {x) satisfy (1.1) and (1.3). Let M > be such that all the 
zeros of all the polynomials det P n (x) are in [-M, M], and let S C C be the set of finite 
cardinality defined in the previous paragraphs. Then for all x G C \ ([— M, M] U S) the 
limiting r x r matrix 

lim P„(x)P- + 1 1 (x) 



exists entrywise and is diagonalizable, with 

lim P n (x)P~K(x)) v fc (x) = z k (x)v k (x) 



k 



1, 



uniformly for k £ {1, . . . , r} and /or x on compact subsets of C \ ([— M, M] U S 1 ). ffere 
we take into account multiplicities as explained in the paragraphs before the statement of 
the theorem. 

Theorem 2.3 is proved in Section 7.3 See also Section [3] for the generalization of 
Theorem 2.3 beyond the matrix Nevai class. 

Since the determinant of a matrix is the product of its eigenvalues, Theorem 2.3 
implies: 



Corollary 2.4. Under the assumptions of Theorem \2.3\ we have 

det P n (x) 



lim 

n->oo det P n+ i(x) 
uniformly on compact subsets ofC\ [-M, M] . 



zi(x) . ..z r (x), 



Note that the convergence in the previous result holds in C \ [— M, M] rather than 
C \ ([—M,M] U S). This is due to Lemma 7.1 See Section [4] below for some further 
corollaries of Theorem 2.3 in terms of the matrix Chebychev polynomials. 



2.3 Limiting zero distribution 

With the above results on ratio asymptotics in place, it is a rather standard routine to 
obtain the limiting zero distribution for the polynomials (P n (x))^ =0 i n matrix Nevai 
class. Recall that the zeros of the matrix polynomial P n (x) are defined as the zeros of 



detP n (x), or equivalently the eigenvalues of the Hermitian matrix J n in (1.2). If these 
zeros are denoted by x\ < X2 < • • • < x rn (taking into account multiplicities), then we 
define the normalized zero counting measure by 

^ rn 

v n = —Y J 5 Xk , (2.12) 
rn ^— ' k 

k=l 

where 5 X is the Dirac measure at the point x. 
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Theorem 2.5. Under the assumptions of Theorem 2.3 the normalized zero counting 
measure u n defined in (2.12) has a weak limit fj,o for n —> oo. The (probability) measure 

tentia, 

xec\r , 



fiQ is supported on the set Tq defined in (2.7) and has logarithmic potential 

J log \x - t\~ l dfi (t) = ^ log \zi (x) . . . z r (x)\ + C, 
for some explicit constant C (actually C = — ^ log | det^4|.j 



(2.13) 



Recall that a measure jiQ on the real line is completely determined from its logarithmic 
potential [25J. 



Theorem 2.5 will be proved in Section |7.4[ with the help of Corollary 2.4 In the 
proof we will obtain a stronger version of ( 2.13[ ), with the absolute value signs in the 
logarithms removed. Moreover, in Section 3] we will extend Theorem 2.5 beyond the 
matrix Nevai class. 

It can be shown that the measure ^o m Theorem 2.5 is absolutely continuous on 
Tq C M with density (see also O [10] ) 



dno(x) 



1 1 

r 2ni 



z' j+ (x) 4_(x) 



z j+ (x) 



dx, 



o- 



(2.14) 



Here the prime denotes the derivation with respect to x, and Zj+(x) and Zj-(x) are the 
boundary values of Zj(x) obtained from the upper and lower part of the complex plane 
respectively. These boundary values exist for all but a finite number of points x 6 IV 
See Section [5] for a comparison to the formulas of |15j . 

2.4 Alternative description of the limiting zero distribution 



First we give an alternative description of the set Tq in (2.7). For the proof see Sec- 
tion EH 



Proposition 2.6. We have 

r = {x G C I 3z e C with f(z,x) = and \z\ = 1} C R. (2.15) 
In words, To is the set of all points x for which f(z, x) = has a root with unit modulus. 
Now let 

X := (x 1 ,x 2 ) C T 

be an open interval disjoint from the set of branch points of the algebraic equation 
f(z,x) = 0. We can then choose a labeling of the roots so that each Zk(x), x G T, is 
the restriction to I of an analytic function defined in an open complex neighborhood 
Note that we do not insist to have the ordering (2.4) anymore. If k is such that 
I Zk (x) I = 1 throughout the interval I then we write 



z k (x) 



o i0 k (x) 



x G 1, 



S 



with 6 k a real valued, differentiable argument function on X. Observe that 



z k (x) 



ie' k (x), 



xGl. 



Moreover 0' k {x) describes how fast z k {x) runs on the unit circle in function of x £ X. 
We can now give an alternative description of the measure no in (2.14) in the following 



proposition which is proved in Section 7.6 



Proposition 2.7. With the above notations we have 



d/ip(x) 
dx 



1 
1 

2irr 



E 

k:\z k (x)\ = l 

E 

k:\z k (x)\=l 



4(*) 



z k (x) 
0' k {x)\ 



x ex. 



(2.16) 
(2.17) 



Moreover, 0' k {x) ^ for any x £ X and for any k with \z k (x)\ = 1. 



3 Generalizations of the matrix Nevai class 
3.1 Slowly varying recurrence coefficients 

In this section we consider a first type of generalization of the matrix Nevai class. We 
will assume that the recurrence coefficients A n and B n depend on an additional variable 
N > 0, as in Dette-Reuther [llj. We write A n ^ and B n ^. For each fixed N > 
we define the matrix- valued polynomials P n> jy(x) generated by the recurrence (compare 
with pTTj )) 

xP njN (x) = A n+ljN P n+ljN (x) + B niN P niN (x) + A* n N P n _ ljN (x) , n > 0, (3.1) 

with again the initial conditions Po.n = Ir an d P—\,n = 0. 
Assume that the following limits exist: 

lim A n , N =: A s , lim B n N =: B s , (3.2) 

n/iV->s n/N^-s 

for each s > 0, where lim n ^^. s means that we let both n, N tend to infinity in such a 
way that the ratio n/N converges to s > 0. We assume that each A s is non-singular. 
We trust that the notation A S ,B S (s > 0) will not lead to confusion with our previous 
usage of A n ,B n (n G Z> ). 

For each s > define the algebraic equation 

= fs(z,x) := det(A*z + B s + A.z' 1 -xl r ). (3.3) 



We again define the roots z k (x,s), k = 1,. . . ,2r to this equation, ordered as in (2.4), 



and we find the corresponding null space vectors v k (x, s) as in (2.2). We also define the 



finite cardinality set S s C C as before, and we define Tq(s) as in (2.7). 



We now formulate the extensions of Theorems |2.3| and |2.5| to the present setting. 
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Theorem 3.1. Fix s > and assume (3.1) and ( |3.2[ ) with the limits A S ,B S depending 
continuously on s > 0. Then for all x G C \ ([— M, M] U 5 S ) £/ie limiting r x r matrix 

lim P n ,7v(a:) J P n " + 1 i i7 v(a;) 
exists entrywise and is diagonalizable, with 

lim P„ i Ar(a;)P~_ 1 1 1 ^(x) ) v fc (x,s) = z fc (x,s)v fc (x,s), fc=l,...,r, 

n/N-ts ' / 

uniformly for x on compact subsets ofC \ ([— M, M] U 5s) . W^e a/so /iaue 

det P n ,jv z) = ^ 

rc/iV->s det P n+ i j Ar(x) 
uniformly on compact subsets ofC\ [-M, M] . 



Theorem 3.2. Under the same assumptions as in Theorem 3.1, the normalized zero 
counting measure of det P n j^{x) for n/N — > s has a weak limit /xo,s; with logarithmic 
potential given by 

f loglx-tr 1 du Jt) = — [ log \zi(x,u) . . . z r (x,u)\ du+C s , x G C\ I J T (u), 
J rs J V,_ 



0<u<t 



for some explicit constant C s 



(3.4) 



In other words, fxo,s is precisely the average (or integral) of the individual limiting 
measures for fixed u, integrated over u G [0, s] . 

Theorems 3.1 and 3.2 are proved in Section 7.7 and an application of these results 
in the context of random band matrices will be given in Section [6j 

3.2 Asymptotically periodic recurrence coefficients 

In this section we consider a second type of generalization of the Nevai class. We will 
assume that the matrices A n and B n have periodic limits with period p G Z>o in the 
sense that 

lim A pn+j =: A®, lim B pn+j =: B® , (3.5) 

n— >oo n— XX) 

for any fixed j = 0, 1, . . . ,p — 1. 

The role that was previously played by zA* + B + z~ l A — xl r will now be played by 
the matrix 



/B^-xl r AW zA®> \ 

A m* bW-xl a& o o 



F(z,x) 














A (p-l) 

^(p- 1 )* ^(p- 1 ) - xl r j 



(3.6) 



prxpr 
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where we abbreviate A&* := (A^)* . We define 

f(z,x) := det F(z, x). 

This can be expanded Laurent polynomial in z: 

r 

f{z,x) = fk{x)z k , 

k=—r 

where now the outermost coefficients f r (x), /_ r (x) are given by 

fr(x) =f r = (-l)^ 1 ^ det(A(°)* . . . A<*-V*), f_ r (x) = f_ r 



(3.7) 



(3.8) 



■ l^y det(A^ ...A^-V) 



We again order the roots Zk{x), k = 1, . . . , 2r of (3.8) as in ( |2.4| ). We define To as in 
(2.7). Proposition 2.1 now takes the following form: 

Proposition 3.3. Tq is a subset of the real line. It is compact and can be written as 
the disjoint union of at most pr intervals. 

As before we denote with v^(x), k = 1, . . . , 2r, the normalized null space vectors of 
the matrix F(zk(x), x). We partition these vectors in blocks as 



Vk(X) 



( v fc,o( x ) 



(3.9) 



2.2 



remains 



where each Vkj(x), j = 0, 1, . . . ,p— 1, is a column vector of length r. Lemma 
true in the present setting. 

Now we can state the analogues of Theorems 2.3 and |2.5| 

Theorem 3.4. For all x G C\ ([— M, M] U S) and any j G {0, 1, . . . ,p— 1}, £/ie limiting 
r x r matrix 



lim P, 



pn+j 



-1 



(X) 



exists entrywise and is diagonalizable, with 

CfiSo P P«+i( x ) P pn+i+p( x )) v fcj( x ) = z k{x)v k ,j (x), (3.10) 
uniformly for k G {1,. . . , r} andforx in compact subsets of C\([— M , M]U S) . Moreover, 

det P n (x) 



lim 



>oo detP n+ p(x) 
uniformly for x in compact subsets of C \ [— M, M] 



zi(aj) . . .Zr(x), 



(3.11) 



Theorem 3.4 is proved in Section 7.8 , where we will establish in fact a stronger variant 
of (3.10). The proof will use some ideas from [3]. Note that the eigenvectors in 
depend on the residue class modulo p, j G {0, 1, . . . ,p — 1}, while the eigenvalues 



7.18 



3.10 



are independent of j. 
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Theorem 3.5. Under the same assumptions as in Theorem \3.4\ the normalized zero 
counting measure of det P n (x) for n — > oo has a weak limit no, supported on To, with 
logarithmic potential given by 

[log\x-t\~ 1 dno{t) = —log\z 1 (x)...z r (x)\+C, xe<C\r , (3.12) 
J pr 

for some explicit constant C (actually C = — log | det(^4 < -°- > . . . j4( p-1 ))|J 

Finally we point out that the results in the present section can be combined with 



those in Section 3.1 That is, one could allow for slowly varying recurrence coefficients 
An,N, -Bn.AT with periodic local limits of period p. The corresponding modifications are 
straightforward and are left to the interested reader. 



4 Formulas for matrix Chebychev polynomials 

Throughout this section we let A be a fixed nonsingular, and B a fixed Hermitian r x r 
matrix. The matrix Chebyshev polynomials of the second kind are defined from the 
recursion (1.1), with the standard initial conditions Po(x) = I r , P-i(x) = 0, and with 
constant recurrence coefficients A n = A and B n = B for all n. The matrix Chebyshev 
polynomials of the first kind are defined in the same way, but now with A\ = ^/2A and 
A n = A for all n > 2. 

Let X and W be the spectral measures for the matrix Chebychev polynomials of the 
first and second kind respectively, normalized in such a way that 



dW = I, 



V ■ 



as in |13t I15j. Here the integrals are taken entrywise. Denote the corresponding Stieltjes 
transforms (or Cauchy transforms) by 

J x-t J x-t 

Theorem 12.31 can be reformulated as 

lim P n {x)P~lAx) = V(x)D(x)V- 1 (x), xeC\ ([-AT, M\ U 5), (4.1) 

n— »oo 

where D{x) is the diagonal matrix with entries Zk(x), k = 1, . . . , r and V{x) is the matrix 
whose columns are the corresponding vectors Vfc(x) in (2.2). 

In turns out that the Stieltjes transforms F\y and Fx can be expressed in terms of 
the matrices D = D(x) and V = V(x) as well. 

Proposition 4.1. We have 

F w (x) = V(x)D(x)V- 1 (x)A- 1 (4.2) 
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and 



F x (x) = [xl r -B-2A*F w {x)A]~ 1 
= V [AVD- 1 - A*VDy l , 



(4.3) 
(4.4) 



for all but finitely many x G C \ Y$. Hence the matrix-valued measures W and X are 
both supported on Tq together with a finite, possibly empty set of mass points on R. 



Proof. By comparing (4.1 ) with Duran's result [131 Thm. 1.1] we get the claimed expres- 
sion (4.2). The formula ( |4.3[ ) follows from the theory of the matrix continued fraction 
expansion, see [2] and also [30] . Formula (4.4) is then a consequence of (4.2)-(4.3) and 
the matrix relation 



(xl r - B)V + A*VD + AVD~ 



0. 



(4.5) 



which is obvious from the definitions of D = D{x) and V = V(x). 



To prove the remaining claims, we note that the right hand side of (4.2) is analytic 



for x G C\(roU5uS') with S the set of branch points of (2.5), so Fyy(x) is also analytic 
there and therefore the measure W has its support in a subset of (Tq US US) PlK. Finally, 



the determinant of the matrix in square brackets in (4.3) is analytic and not identically 



zero for x G C \ (ro U5U5) so it has only finitely many zeros there. This yields the 
claim about the support of the measure X. □ 

The above descriptions considerably simplify if A is Hermitian. In that case the 



algebraic equation (2.5) becomes 



= f(z, x) = det(Au; + B - xl r 



where 



w 



z + z 



(4.6) 



(4.7) 



As in Section 2.4, for any open interval I := (^1,^2) C Tq disjoint from the set of 



branch points of the algebraic equation (2.5), we can choose an ordering of the roots 



Zk(x), x G X, so that each Zk(x) is the restriction to I of an analytic function defined on 



an open complex neighborhood ft D X. Thereby we drop the ordering constraint (2.4) 

1 r 



By (4.6)-(4.7) we may assume that 



Z2r-k(x) = Z k {x) \ 

for all x G fl D I and we write 

w k (x) = z k {x) + z k (x)~ 
If w k {x) G (—2, 2) for all x G X then write 
w k (x) = 2cos9 k (x), 



k 



1. 



< 9 k (x) < tv, 



(4.8) 



(4.9) 



(4.10) 



with 9 k a real-valued, differentiable argument function on X. Denote with V{x) the 
matrix formed by the normalized null space vectors v^(x) for the roots w k (x), x G X. 
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Proposition 4.2. Assume that A is Hermitian. Then with the above notations, the 
density of the absolutely continuous part of the measures X and W is given by 



dX(x) 
dx 

dW(x) 
dx 



V(x)Ax(x)V- l (x)A~\ x G r , 



V(x)A w (x)V- 1 (x)A-\ 



xer, 



o, 



where 



1 



Ax(z) = -diag 



L«, fc (x)e(-2,2) 



sign w' k (x) 



k=l 



v y/A ~ W k {x) 2 

&w(x) = ^ dia § ( 1 «'feWe(-2,2)V /4 - w k {x) 2 s\gnw' k {x) 



r 

k=l 



and where the characteristic function l Wke r_ 2j 2) takes the value 1 if w k 6 (—2,2) and 
zero otherwise. 



Proposition 4.2 is proved in Section 7.9 If A is positive definite then the factors 
signu^(x) in the above formulas can be removed. This follows from (5.2) below. Then 
one can show that the above formulas correspond to those in |13} 115] . 



5 The results of Duran-Lopez-Saff revisited. 



In this section we show how Theorem [2^5] on the limiting zero distribution of P n (x) in the 
matrix Nevai class, relates to the formulas of Duran-Lopez-Saff [15] for the case where 
the matrix A is positive definite or Hermitian. 

Throughout this section we write the algebraic equation f(z, x) = as in (4.6)— (4.7). 
First we will assume that A is positive definite. Then A 1 / 2 exists and we can replace 
the algebraic equation (4.6) by 

= det(W r + A- l ' 2 BA- 1 ' 2 - xA' 1 ). 

Hence the roots w are the eigenvalues of the matrix 

xA -i _ A-V 2 BA-V 2 . 

If x € R then this matrix is Hermitian and we denote its spectral decomposition by 

xA- 1 - A-WBA-V 2 = U(x)D w {x)U- 1 (x), (5.1) 

where 

D w (x) = diag(wi(x), ... ,w r (x)) 

is the diagonal matrix containing the eigenvalues, and U(x) is the corresponding eigen- 
vector matrix. We can assume that U(x) is unitary, i.e. U^ 1 (x) = U*(x). 

As in the previous section, we fix an open interval X := (x\,X2) C Tq disjoint from 
the set of branch points of the algebraic equation (2.5), and we choose an ordering of 
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the roots z k (x), x £ I, so that each z k (x) is the restriction to X of an analytic function 
defined on an open complex neighborhood 0, D I. The same then holds for w k (x) in 

Lemma 5.1. If A is positive definite, then with the above notations we have 

w' k (x) = (U- l (x)A- 1 U(x)) k>k >0, k = l,...,p, (5.2) 
where we use the notation M k) k to denote the (k, k) entry of a matrix M . 



Proof. We take the derivative of (5.1) with respect to x. This yields (we suppress the 
dependence on x for simplicity) 

A' 1 = XJ'D W \J- X + UD'JJ- 1 + IID^U- 1 )', 

or 

U- x A- l U = \J- X XJ'D W + D' w + D w {U' lv 
Invoking the fact that (U^ 1 )' = —U~ l U'U~ l , this becomes 

XJ- X A- X XJ = D' w + [U- l U\D w ] , 



U. 



where the square brackets denote the commutator. The equality in (5.2) then follows 
on taking the (k, k) diagonal entry of this matrix relation, and noting that the diagonal 
entries of the commutator \U~ l U' , D w ] are all zero since D w is diagonal. Finally, the 
inequality in (5.2) follows because A is positive definite and U is unitary, U^ 1 = U*. □ 



From (4.9) we have 



40*0 = ±i w 'k( x ) 

Zk(x) yjA-w k {xf 



Thus the density of the limiting zero distribution fiQ in (2.16) becomes 



d^ojx) 
dx 



w'Jx) 



1 r 



We(-2,2), 



(5.3) 



(5.4) 



where we used that w' k (x) > 0, and where again the characteristic function 1^ g(_2,2) 
takes the value 1 if Wk € (—2,2) and zero otherwise. Note that the factor 2 in the 
denominator of (2.16) is canceled since for any w k 6 (—2,2) there are two solutions z k 
to (4.9), one leading to the plus and the other to the minus sign in (5.3). 



Inserting (5.2) in ( |5.4[ ) we get 

d/j, (x) 



1 



dx 



irr 



E 

k=l 



(U~ 1 (x)A- 1 U(x)) ktk _ 



w k (x) 



-Tr (U- 1 (x)A- 1 U(x)A x (x)) 
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where Tr(C) denotes the trace of the matrix C and Ax(x) is the diagonal matrix 



Ax(*) = -diag( W( ~ 2 ' 2) 
vr \y/4:-w k (x) 2 



k=l 



Since the trace of a matrix product is invariant under cyclic permutations, we find 

= -Tr(A- 1 / 2 U(x)A x (x)U- 1 (x)A-y 2 ) . 
dx r \ J 

This is the result of Duran-Lopez-Saff for the case where A is positive definite |15j . 

Next suppose that A is Hermitian but not necessarily positive definite. Then we can 
basically repeat the above procedure. Rather than taking the square root A 1 / 2 , we now 



write (4.6) in the form 



= det(wl r + A~ X B - xA' 1 ). 

Hence the roots w are the eigenvalues of the matrix xA -1 — A~ 1 B. Supposing that this 
matrix is diagonalizable, we write 

xA- l -A~ 1 B = V{x)D w {x)V~ 1 {x), (5.5) 

with D w = diag(wi, . . . , w r ) but with V not necessarily unitary. The notation V is 



compatible with the one in the previous section, by virtue of (4.5). We then again have 



the equality in (5.2) (with U replaced by V) although the positivity w' k {x) > may be 
violated. Similarly to the above paragraphs, we find 



d/j, (x) 



-Tr (V(x)Ax(x)V~ 1 (x)A- 1 ) , 



dx r 

with Ax defined in Proposition |4.2| The latter proposition then implies that 

d^o(x) = l_ Tr f dX(x) 
dx r \ dx 

This is the result of Duran-Lopez-Saff for the case where A is Hermitian [15 



6 An application to random matrices 

In a fundamental paper Dumitriu-Edelman [12] introduced a tridiagonal random matrix 



/ 



G (i) 



Ni 



\ 



N 3 



n-l 



V2 



Xr 



16 



where Ni,..., N n are independent identically standard normal distributed random vari- 
ables and Xip, . . . , are independent random variables also independent of Ni, . . . , 
N n , such that Xf g is a chi-square distribution with jfi degrees of freedom. They showed 



that the density of the eigenvalues Ai < • • • < A n of the matrix G r , 
called beta ensemble 

n ^2 

^ni^- A ii^ ex p(-E^ 

i<j j=l 



is given by the so 



where cp > is an appropriate normalizing constant (see [16J or [22j among many 
others). It is well known that the empirical eigenvalue distribution of ~^G n converges 
weakly (almost surely) to Wigner's semi-circle law. In the following discussion we will 



use the results of Section 3.1 to derive a corresponding result for a (2r + l)-band matrix 



of a similar structure. To be precise consider the matrix 



/ \/2 Ni A'(„_ 1 ) 71 



X( n —2)j 2 X( n — 2)7! 



V2 



X, 



(n — ^J'Yt- 



X, 



(n-r)7 r _! 



Xi, 



-l)7r 



A', 



(n — 3)72 



X, 



(n — r — l)^ r 



x 2jl 



X2J1 



X 12 

Xyx 

V2N n ) 



(r) 

where n = mr, all random variables in the matrix G n are independent and Nj is 
standard normal distributed (j = 1, . . . , n), while for j = 1, . . . , n — 1, k = 1, . . . , r 
has a chi-square distribution with j^ k degrees of freedom (71, . . . , j r > 0). It now 
follows by similar arguments as in [11] that the empirical distribution of the eigenvalues 



A 



(n,r) 



< 



< Xn of the matrix 



(r) 

Gn has the same asymptotic properties as the 



limiting distribution of the roots of the orthogonal matrix polynomials R mtn (x) defined 
by (i?_i >n (x) = 0, R 0) n(x) = I r ) 

xR k>n {x) = A k+1>n R k+ x >n (x) + B k>n R k>n (x) + A* k ^ n R k _ ltn (x); k>0, 
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where the r x r matrices and Bi^ n are given by 

/ - l)r + l)jr 

V((i-l)r + 2) 7r _ 1 ^((i-l) r + 2) 7r 



.4 



t.ri 



5; 



1 



\ 





V 



\J {i r - 1 )72 



\J{ir- l)7r-i \/(zr - l) 7r 



( x/(ir + l) 7 i + 1)72 

ypr + 1)71 ^{ir + 2)71 



y /gr + l)7r-2 

V \/(ir + l) 7r -i 



V((t + l)r - 2)71 



x /((» + lV-2) 72 v^i + l)r - l) 7l 
Now it is easy to see that for any u £ (0, 1) 

lim B ijTl = B(u), lim A i<n = A(u), 



yf{ir + i) 7r _ 1 \ 
v/pr + 2) 7r _ 2 



+ l)r - 2) 7l 




where 



\/7r-l jl~r 



72 



/ 



y /lr-2 
\ y/lr-1 



\ 




V7r-1 v 7 ^ 

y/lr-2 y/7r-l jlr / 



71 V72 
71 



v 7 ^ ^/tI 

V 7 ^ v 7 ^! / 



6 



(6.1) 



(6.2) 



and Theorem |3.2| shows that the normalized counting measure of the roots of the ma- 
trix orthogonal polynomial R m ^ n {x) has a weak limit /u 0i i/ r defined by its logarithmic 
potential, that is 



log|x-t| dno,i/ r {t) 



l/r 



log \ zi(x, u), . . . , z r (x, u)\du + c(r) 
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where z\(x,u), . . . , z r (x, u) are the roots of the equation 

det(A*(u)z + B(u) + A^z" 1 - xl r ) 







(6.3) 



corresponding to the smallest moduli. Observing the structure of the matrices in (6.1) 
and (6.2) it follows that 



where zi(x) , . . . , z r (x) are the roots of the equation (6.3) for u = 1 (with smallest 
modulus). Similar arguments as given in the proof of Proposition 2.7 now show that the 
measure /xo,i/r is absolutely continuous with density defined by 



d/4),l/r0) 



d.r 



1/r 1 
2^ 



1 

2^ 



E 

fe:|*fc(a:,«)[=l 
1/r j 



Ox 



z k (x,u) 



du 



E 

fc:|**(a!/-v/«)|=l 



z k (x,u 



Zk(x/Vu) 



du. 



(6.4) 
(6.5) 



By the same reasoning as in [TT] we therefore obtain the following result (5 X denotes the 
Dirac measure). 



Theorem 6.1. If X^' r ^ < ••• < An"'''' denote the eigenvalues of the random matrix 
Gn I \f n ~ with 71, . . . ,7 r > 0, then the empirical eigenvalue distribution 

1 " 



, (n,r) 



converges weakly (almost surely) to the measure /^o,i/r defined in (6.4)- 

We conclude this section with a brief example illustrating Theorem 6.1 in the case 
r = 2. In the left part of Figure [l] we show the simulated eigenvalue distribution of the 

( r ) , , 

matrix Gn /yn in the case n = 5000 71 = 72 = 1, while the corresponding limiting 
density is shown in the right part of the figure. Similar results in the case r = 2, 71 = 1 
and 72 = 5 are depicted in Figure [2J Note that the derivatives z' k (x) can be evaluated 
numerically using the formula (implicit function theorem) 



H( z k(x),x) 



7 Proofs 



7.1 Proof of Proposition [27T 



First we consider the behavior of the roots z^ix) for x — > 00. It is easy to check that in 
this limit, 

z k (x)^0, k = l,...,r, 

\zk(%)\ ->■ 00, k = r + 1, . . . ,2r, 



(7.1) 
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Figure 1: Simulated and limiting spectral density of the random block matrix G n / 1 y/n 
in the case r = 2, 71 = 1, 72 = 1. In the simulation the eigenvalue distribution of a 
5000 x 5000 matrix was calculated (i.e. m = n/r = 2500). 




E 6 



Figure 2: Simulated and limiting spectral density of the random block matrix Gn / 1 \fn 
in the case r = 2, 71 = 1, 72 = 5. In the simulation the eigenvalue distribution of a 
5000 x 5000 matrix was calculated (i.e. m = n/r = 2500). 
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see e.g. [10] . In particular we have that 1^(^)1 < |<Sr+i(a;)| if \x\ is large enough, showing 
the compactness of IV 

Next we ask the question: for which ieC can we have a root z k = Zk(x) such that 
\zk\ = 1? In this case, (2.5) becomes 

= f(z k , x) = det(z k A* +B + zj,A- xl r ), 

so x is an eigenvalue of the Hermitian matrix z^A* + B + z^A. In particular, it follows 
that 

|z fc (jc)| = 1 => xeK. (7.2) 



Now for x — > oo, (7.1) implies that we have precisely r roots Zk(x) with |zfc(x)| < 1. By 
(7.2) and continuity this must then hold for all x G C \ M: 

\z k (x)\ < 1, k = l,...,r, 



x e c\: 



> 1, fc = r + 1, 



,2r. 



In particular we see that |z r (a;)| < 1 < |z r +i(a;)| for i£C\l, implying that To C R. 

Finally, the claim that To C IR is the disjoint union of at most r intervals follows 
from HE [2S| . □ 



7.2 Proof of Lemma EO 

The proof will use ideas from Duran [131 Proof of Lemma 2.2]. We invoke a general result 
known as the square-free factorization for multivariate polynomials. In our context it 
implies that there exists a factorization of the bivariate polynomial z r f(z, x) of the form 

K 

z r f(z,x) = l\g k (z,xr\ (7.3) 

k=l 

for certain K G Z>o, multiplicities mi, . . . ,itik £ ^>o an d non-constant bivariate poly- 
nomials gi(z, x), . . . , gx(z, x), in such a way that 

K 

g{z,x) := Y\_gk(z,x) (7.4) 

k=l 

is square-free, i.e., for all but finitely many x £ C, the roots z of g(z,x) = are all 
distinct, and vice versa with the roles of x and z reversed. 

The existence of the square-free factorization of the previous paragraph, can be 
obtained from a repeated use of Euclid's algorithm. For example, if z r f(z,x) = has 
a multiple root z = z(x) for all x £ C, then we apply Euclid's algorithm with input 
polynomials z r f(z,x) and ^ (z r f(z,x)), viewed as polynomials in z with coefficients in 
C[x}. This gives us the greatest common divisor of these two polynomials and yields a 
factorization 

z r f(z,x) = gx(z,x)g 2 (z,x), 
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for two bivariate polynomials gi,g2 which depend both nontrivially on z. Note that the 
factorization can be taken fraction free, i.e., with g±,g2 being polynomials in x rather 
than rational functions. If one of the factors g\ or g<± has a multiple root z = z{x) for all 
x € C, then we repeat the above procedure. If g\ and 52 have a common root z = z(x) 
for all x £ C, then we apply Euclid's algorithm with input polynomials g\ and (72, viewed 
again as polynomials in z with coefficients in C [x] . Repeating this procedure sufficiently 
many times yields the square-free factorization in the required form. 



Note that the factors g^{z,x) in (7.3) all depend non-trivially on both z and x. For 



if gk(z,x) would be a polynomial in z alone (say), then there would exist z £ C such 



that f(z,x) = for all x G C, which is easily seen to contradict with (2.5). 
From the above paragraphs we easily get the symmetry relation 

mi(z, x) = m2(z,x), 



for all but finitely many z, x £ C. This proves one part of Lemma 2.2 

It remains to show that m2(z, x) = d(z,x) for all but finitely many z,x € C. From 
the definitions, this is equivalent to showing that the matrix zA* + B + z~ l A is diag- 
onalizable for all but finitely many z £ C. This is certainly true if \z\ = 1 since then 
zA* + B + z~ 1 A is Hermitian so in particular diagonalizable. The claim then follows in 
exactly the same way as in [131 Proof of Lemma 2.2]. □ 



7.3 Proof of Theorem [231 



Before going to the proof of Theorem 2.3 let us recall the following result of Dette- 
Reuther As mentioned before, this result uses in an essential way the Hermitian 



structure of (1.2). 

Lemma 7.1. (See fllf:) If all roots of the matrix orthogonal polynomials P n {x) are 
located in the interval [—M,M], then the inequality 

W*Pn(z)P-^(z)A-l lV \ < 



distO, [-M, M\) 

holds for all complex numbers z and for all column vectors v with unit Euclidean norm 
llvll = 1. 



We will need the following variant of Lemma 7.1 



Corollary 7.2. In the matrix Nevai class (1.3) we can find M > as in Lemma 7.1 
that 

8\\A\\ 



so 



|w*F»WP„-iWv|< 5ist(2|[ _ M>M]) 



(7.5) 



for all n sufficiently large, all column vectors v,w with ||v|| = ||w|| = 1 and all complex 
numbers z. Here \\A\\ denotes the operator norm (also known as 2-norm, or maximal 
singular value) of the matrix A. 
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Proof. For fixed n and fixed z € C, define the sesquilinear form 

(v,w) A := w*P n (z)p-^ 1 (z)A-l 1 v, 

which is linear in its first and antilinear in its second argument. Also define ||v||^ := 
(v, v)a- (This 'norm' is not necessarily positive!) The polar identity asserts that 



(v,w) A = - (||v + w|| A 



|v — w|| A + i||v + iw" z 



A -i\\v-iw\\ 2 A ) 



Combining this with Lemma 7.1 we get 

|w*P n ( 2 )p- + 1 1 (. ; )^ 1 v| =: |(v,w) A | < 



dist(z, [-M, M\) 



for all pairs of vectors v,w with ||v|| = ||w|| = 1 and for all complex numbers z. If we 
now take n sufficiently large so that ||j4 n+ i|| < 2||A|| (recall ( |1.3[ )), we get the desired 
inequality (7.5). □ 

Proof of Theorem \2.3[ We will use a normal family argument. Fix k £ {1, ... ,r} and 
zq € C \ [—M,M]. By (7.5) we have that P n (z)P~+ 1 (z)v k (z) is uniformly bounded 
entrywise in a neighborhood of z = oo. By Montel's theorem we can take a subsequence 
of indices (nj)°^ so that the limit lim^oo P ni (z)P~ +1 (z)v k (z) exists uniformly in this 
neighborhood. We will prove by induction on I = 0, 1, 2, . . . that 



lim P ni (x)P. ' i Jx)v k (x) = z k (x)v k (x) + 0(x l ), 



x — > oo. 



(7.6) 



For I = (or I = 1) this follows from (7.5) and the fact that z k {x) = 0(x~ 1 ) as x — > oo. 

Now assume that (7.6) is satisfied for a certain value of I, and for any sequence (n«)j 
for which the limit exists. Fixing such a sequence (rtj)j, we will prove that (7.6) holds 
with / replaced by I + 2. By moving to a subsequence of {rii)i if necessary, we may 
assume without loss of generality that the limiting matrices lim^oo P ni (x)P ~ +1 ( x) and 
linxj-^jx) P-ns— 

i(x)P n 1 (x) both exist. The induction hypothesis asserts that ( |7.6[ ) holds 
for any sequence (nj) for which the limit exists, so in particular 



lim P ni _i(x)P 1 (x)v fc (x) = z k {x)v k (x) + 0(x l ), 



X — > CO. 



(7.7) 



Now let us write the three-term recurrence (1.1) in the form 

A^Pn^P-^x) + {B n . - Xl r ) + A ni+l P ni+l (x)P^(x) = 0. 

Multiplying on the right with v k (x), taking the limit i — > oo and using the facts that 
lim n A n = A and lim n B n = B we get 

A* lim (P ni _ 1 (x)P" 1 (x)) v k {x) + (B -xl r )v k (x) + A lim (P rH+l {x)P^{x)) v k {x) = 0. 

With the help of ((7^| and ((2^2} this implies 

A lim (P n . + i(x)p- 1 (x)) v k (x) = Az^ix^kix) + 0(x~ l ), x -> oo. 
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Multiplying this relation on the left with z k (x) lim n _ s>00 [P Tli (x)P n . 1 +1 (x)) A 1 and using 
(7.5) and the fact that z k (x) = 0(x^ 1 ) for x — > oo, we find 

lim P n .{x)P~\ 1 {x)v k {x) = z k (x)v k (x) + 0(x^ 2 ), x -)■ oo, 



showing that the induction hypothesis (7.6) holds with I replaced with 1 + 2. This proves 
the induction step. □ 

7.4 Proof of Theorem 12.51 

We write the telescoping product 

detP„(x) = (detP n (x)p-\{x)) ...(detP 2 (x)P{ 1 (x)) (det P 1 (x)Pq 1 (x)) detP (x). 
Taking logarithms and dividing by rn we get 



1 log det P n (x) = — ( ( V log det P k {x)P k -\ (x) ) + log det P (a 
m rn \ \ ^ / 



(Here we use the logarithm as a complex multi-valued function.) Taking the limit n — > oo 
and using the ratio asymptotics in Corollary 2.4, we obtain 

- lim — log( det PJx)) = -log(z 1 (x) ...Zr(x)), x G C \ \—M, Ml. 
n— >oo rn r 

Now we take the real part of both sides of this equation. Then the left hand side becomes 
precisely the logarithmic potential of fiQ, up to an additive constant C. So we obtain 
(2.13); the constant C can be determined by calculating the asymptotics for x — > oo. □ 



7.5 Proof of Proposition 2.6 



From Proposition 2.1 and its proof we know that both the left and right hand sides of 
the equation in (2.15) are subsets of the real axis. Now for i£l one has the symmetry 
relation 



f(z, x) := det(A*z + B + Az~ x - xl r ) = det(Az + B + A*z~ L - xl r ) = f(z~ 



(7.8) 



where the bar denotes the complex conjugation and we used that detM = det(M*) for 
any square matrix M. This implies that for each solution z = z k {x) of the equation 
f(z,x) = 0, the complex conjugated inverse z = z^, l {x) is a solution as well, with the 



same multiplicity. So with the ordering \2A\ we have that |#jfe(a?)| ■ \z2 r -k{x)\ = 1 for 
any k = 1, . . . ,r. In particular we have that |z r (a;)| = |z r+ i(x)| if and only if |2 r (x)| = 1. 
This implies Proposition 2.6 □ 
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7.6 Proof of Proposition 



2.7 



We use the notations in Section 12.41 We fix x 6 X C Tn and define the sets 



S+(a:) = {z k (x) | |^fe(y)| < 1 for all y € O n C+}, 
S-(x) = {z k (x) \ \z k (y)\ < 1 for all y efinC-}. 



(7.9) 
(7.10) 



So S + (x) (or £_(„)) contains all the roots z k (x) for which |zfc(z/)| < 1 for y in the upper 
half plane (or lower half plane respectively) close to x. 

Let Zk(x) be a root of modulus strictly less than 1. By continuity this root belongs 
to both sets S+(x) and S-(x), with the same multiplicity, and hence the contributions 
from the +- and — terms in ( |2.14 ) corresponding to this root z k (x) cancel out. 

Next let z k (x) be a root of modulus 1. Assume again that z k (y) = e i0k ^ with k (y) 
real and differentiable for y G Z C M. Suppose that k (x) > 0. Then the Cauchy- 
Riemann equations applied to log z k (y) imply that | £&(_)! < 1 for y in the upper half 
plane close to x, and |£fc(y)| > 1 for y in the lower half plane close to x. So z k {x) lies in 
the set S + {x) in (7.9) but not in S-(x). Similarly if 0' k (x) < then z k (x) lies in the set 
S-(x) but not in S+(x). In both cases, the contribution from z k (x) in the right hand 



side of (2.14) has a positive sign and so we obtain the desired equality (2.16). 

Finally, the claim that 9'{x) ^ for any x G X C Tq follows since, if this fails, 
then general considerations (e.g. in |26t Proof of Theorem ll.l.l(v)]) would imply that 
Tq <f_ R, which is a contradiction. □ 



7.7 Proof of Theorem I37T1 



The proof of Theorem 2.3 and 2.5 can be easily extended to prove Theorem 3.1 The 



difference is that the limits limn._j.oo should be replaced by local limits of the form 
lim n ./jv^_ The details are straightforward and left to the reader (for similar reasonings 
see also [3 HH DH EI] , among others.) 

7.8 Proof of Theorem 13^1 



The proof of Theorem 3.4 will follow the same scheme as the proof of Theorem |2.3[ but 
it will be more complicated due to the higher periodicity. To deal with the periodicity 
we will use some ideas from [3|. It is convenient to substitute z = y p and work with the 
transformed matrix 



G(y,x) := diag(l,y, . . . ,y p 1 )F(y p , x) diag(l, y \ 

/£?(°) - xl r y~ l A^ 
y AW* - xl r y- l A^ 

yA^> 



V y- x A<® 














y -l A (p-l) 

yAte-V* B(P~V - xl r J 



(7.11) 



(7.12) 



prxpr 
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Consistently with the substitution z = y p , we put y k (x) = z¥ p (x), k = 1, . . . , 2r, for an 
arbitrary but fixed choice of the pth root. The ordering (2.4) implies that 

< \yi{x)\ < \y 2 (x)\ <...< \y r (x)\ < \y r +i{x)\ <...< \y 2 r{x)\. (7.13) 

Note that each y = yk{x) is a root of the algebraic equation 

det G(y, x) = det F(y p , x) = 0. 



From (7.12) it is then easy to check that (see e.g. |10j ) 

" 1 , k = 1 r, 



Vk(x) oc 



x K = i, . . . , r, x — > oo, 

x, ft = r + 1, . . . , 2r, x — ^ oo, 



(7.14) 



where the oc symbol means that the ratio of the left and right hand sides is bounded 
both from below and above in absolute value when x — > oo. 
Denote with w k (x) a normalized null space vector such that 

G{y k (x),x)w k (x) = 0. (7.15) 

If there are roots yk(x) with higher multiplicities then we pick the vectors Wfc(x) as 
explained in Section 2.2 We again partition w/%(x) in blocks as 

/ w fcj0 (x) 



w fc (x) 



(7.16) 



\w fci p_i(x) / 



where each Wkj(x), j = 0, 1, ... ,p — 1, is a column vector of length r. Assuming the 
normalization ||wfc(x)|| = 1 then we have that 



lim ||w fe j(x)|| = C k ,j > 0, 



j = 0,l,...,p-l. 



(7.17) 



This follows from (7.15)-(7.16), (7.14) and by inspecting the dominant terms for x — > oo 
in the matrix (7.12). 

Theorem 3.4 will be a consequence of the following stronger statement: 

\J^ p pn+j(x)P^ + j +1 (x) \ w kJ+1 (x) = y k (x)w kJ (x), x -)• oo, (7.18) 

uniformly for x in compact subsets of C \ ([— M, M] US), for all k G {1, . . . , r} and 
for all residue classes j E {0, 1, . . . ,p — 1} modulo p. (We identify w kjP (x) = ~w k! o(x).) 
Indeed, Theorem 3.4 immediately follows by iterating (7.18) p times and using that 
y k (x) = z k{x). 

The rest of the proof is devoted to establishing (7.18). We will show by induction on 
I > that 

P P n *+J'( x ) P p^+j+i( x )) ™k,j+i(x) = y k (x)w kd (x)(l + 0(x~ 1 )), x->oo, (7.19) 
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for any ft € {1, ... ,r} and j E {0, 1, . . . ,p — 1}, and for any increasing sequence (nj)°^ 
for which the limit in the left hand side exists. 

Assume that the induction hypothesis ( 7.19| ) holds for a certain value of I > 0. We 
will show that it also holds for I + 2. Let (Vii)^n be an increasing sequence for which 



the limit in the left hand side of (7.19) exists. We can assume without loss of generality 



that j = p — 1; a similar argument will work for the other values of j £ {0, 1, ... ,p — 1}. 
Now from the three-term recursion we obtain 



(A 



V'H 



At 



Pprii (x) 



A 



prii+p-1 



A* R 

Siprii+p-l ^prii+p-l 



xl r Ap ni ^.p J 



Pprii+p-l{% 
\ Ppni+p{x) J 



0. 



Applying a diagonal multiplication with appropriate powers of y := yk(x) we get 



(7.20) 



(yA 



pit, 



B 



pit, 



xl r y Aprn+i 



yA 



prii+l 



y-'A 



yA, 



prii+p-l 



prii+p—1 



B 



pni+p-1 



xl r y Ap ni +p J 



( y-PPpn^ix) \ 
y-^P P n t (x) 

y^Ppm+p-iix) 

Pprii+p-l{x) 
\ yPpn l+ p( X ) J 

(7.21) 



Let us focus on the rightmost matrix in the left hand side of (7.21). Multiplying on the 
right with P^ +p _ 1 (x)w fcp _ 1 (x) it becomes 



/ y P Ppn. l -\Ppn l +p-l^k,p-l \ 

y~ 1 Ppn l +p-2 Ppn, +p- 1 Wfc,p_ 1 
\ yPpni+pPp rH+ p-iWk,p-l / 



(7.22) 



(Here we skip the x-dependence for notational simplicity.) By moving to a subsequence 



of (n«)^Q if necessary and using compactness, we may assume that each block of (7.22) 



has a limit for i — > oo. Repeated application of the induction hypothesis (7.19) then 



implies that the limit of (7.22) for i — > oo behaves as 



fw kjP ^(l + 0(x~ l ))\ 
w kfi (l + 0(x- 1 )) 

w fciP _ 2 (l + 0(x- i )) 

W fe>p _l 
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for x — > oo, where 



<p(x) := ( lim Ppn l + P {x)P p ^ +p _ 1 (x) ) y k (x)-w k:P ^{x). 



(7.23) 



Multiplying (7.21) on the right with P p n i+p ^i{x)wk,p-i{ x ) and taking the limit i 



oo, we get from the above observations that 
f yA (°> B^-xl r y~ l A^ 

yAW* 



y-l A {V-l) 



\ 



yAiP- 1 )* £(p-!) - xl r y~ l A^ j 



f\Vk,p-l\ 
Wfen 



fO{x- l+1 )\ 
0{x- l+1 ) 

yoix- 1 - 1 )/ 



(7.24) 



for x —7- oo, where we used that y = y k (x) oc x 1 for x — >■ oo. Taking the last block row 
of this equation yields 



Hi: 



(x)A^ 1 >w kjP - 2 + (B^-V - xJ r )w fciP _i + y k \x)A^ v (x) = 0{ X - 1 - 1 ). (7.25) 



On the other hand, by ( |7.15| ), ( |7.16 ) and (7.12) (evaluated for the last block row) we 
have that 

y k \x)A^w k>0 + y k {x)A^- l >vt Kp „ 2 + (B^ - xJ r )w fciP _! = 0. 



Subtracting this from (17.25 ) we get 



A® f Um P P n i+P (x)P pn ; +p _ 1 (x)) w^ifr) - ^(x^w^x) = O^ 1 ), 
on account of (7.23). The factor ^4^°^ can be skipped from this equation. Then multi- 



plying on the left with y k (x) x (lim;-^ P prH+p ^ 1 (x)P pr ^ +p (x)) we get 

li^P pn . +p ^i(x)P p ^ i+p (x) \ w M (x) - y k (x)w k j>-i{x) = <3(x~^ 3 ), 
or equivalently 

(j^Pmi+P-i( x ) P mi+ P ( x )) w fc,o( x ) = yk(x)w k , p -i(x)(l + 0(x^~ 2 )). 



We conclude that (7.19) holds with I replaced by 1 + 2. This proves the induction step. □ 
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7.9 Proof of Proposition |4.2 



Throughout the proof we will use the notations of Section [4j Recall that the Hermitian 
symmetry A = A* implies the roots Zk(x) to appear in pairs {zk(x), ^(x) -1 }. Both 
Zk(x) and Zfc(x) -1 correspond to the same value of Wk(x) = Zf~(x) + z k {x)~ l in (4.9) and 
therefore to the same null space vector Vfc(x). 

Now let x G R. For any w k (x) = 2cos9 k (x) G (-2,2), with 6 G (0, ir), k = 1, . . . , r, 
we have a pair of roots 2^ (x) = e ldk ^ and Zfc 2 (x) = e~' l9k ^ x \ with ^2} = {k, 2r — k}. 
Suppose that vo' k (x) > 0. Then the Cauchy-Riemann equations show that z kl [x) lies in 



the set S—(x) in (7.10) but not in S+(x), and vice versa for the root z k2 (x). The reverse 



situation occurs if w' k (x) < 0. 

Fix i£l and assume the labeling of roots is such that 



max{|zi(x)|, . . . , \z K (x)\] < 1, 



\z K +i{x) 



\z r (x)\ = 1, 



with K G {0, ...,r}. Taking into account the above observations, we find from (4.2) 
that 

lim Fw(x + ei) 
= V{x) ( lim + D(x + ei)^ V~\x)A~ x 

= V(x) diag(z!(x), . . . , Z K ( X ), e -^+lWsign^ +1 (x) ; _ ^ e ~ie r (x) S i g n w' r (x) ( x ) A ~l . 



Similarly 



lim i*V(x — ei) 

e->-0+ 



-1 



= V(x) I lim + D(x - ei) 1 V _1 (x)j4. 
= V(x) diag(^(x), . . . , zjr(x), e*+^) si s n <+iW, . . . , ^M^ig^O*))]/-!^-!. 
Using the Stieltjes inversion principle 
dW(x) 1 



dx 



^ . lim (Fw(x — ei) — Fw(x + ei)) , 
2m e-K)+ 



the desired formula for dPF/dx now follows from a straightforward calculation. The for- 

ue 

□ 



mula for dX/dx similarly follows from (4.4), taking into account the simplifications due 
to A = A*. 
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